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Abstract 

We have built a 12DOF, passive-compliant legged, 
tailed biped actuated by four brushless DC motors. We 
anticipate that this machine will achieve varied modes of 
quasistatic and dynamic balance, enabling a broad range 
of locomotion tasks including sitting, standing, walking, 
hopping, running, turning, leaping, and more. Achiev¬ 
ing this diversity of behavior with a single under-actuated 
body, requires a correspondingly diverse array of con¬ 
trollers, motivating our interest in compositional tech¬ 
niques that promote mixing and reuse of a relatively few 
base constituents to achieve a combinatorially growing ar¬ 
ray of available choices. Here we report on the develop¬ 
ment of one important example of such a behavioral pro¬ 
gramming method, the construction of a novel monopedal 
sagittal plane hopping gait through parallel composition 
of four decoupled IDOF base controllers. 

For this example behavior, the legs are locked in phase 
and the body is fastened to a boom to restrict motion to 
the sagittal plane. The platform’s locomotion is powered 
by the hip motor that adjusts leg touchdown angle in flight 
and balance in stance, along with a tail motor that adjusts 
body shape in flight and drives energy into the passive 
leg shank spring during stance. The motor control sig¬ 
nals arise from the application in parallel of four simple, 
completely decoupled IDOF feedback laws that provably 
stabilize in isolation four corresponding IDOF abstract 
reference plants. Each of these abstract IDOF closed 
loop dynamics represents some simple but crucial specific 
component of the locomotion task at hand. We present a 
partial proof of correctness for this parallel composition 
of template reference systems along with data from the 
physical platform suggesting these templates are anchored 

*‘Electrical and Systems Engineering, University of Pennsylvania, 
Philadelphia, PA, USA. {avik, kod}@seas. upenn. edu. 

^This work was supported in part hy the ARL/GDRS RCTA 
project. Coop. Agreement #W91 lNP-1020016 and in part by NSF 
grant #1028237. 


as evidenced by the correspondence of their characteristic 
motions with a suitably transformed image of traces from 
the physical platform. 

1 Introduction 

The control of power-autonomous, dynamic legged 
robots that have a high number of degrees of freedom 
(DOF) is made difficult by a number of factors including 
(a) under-actuation necessitated by power-density con¬ 
straints, (b) the existence of significant inertial coupling 
and Coriolis forces that are hard or impossible to can¬ 
cel, (c) variable ground affordance, (d) often hard-to- 
measure and necessarily rapid hybrid transitions. In the 
face of these challenges, some popular methods of con¬ 
troller design, such as hybrid zero dynamics [2]— which 
are “exact” in their domain of applicability but require 
extremely accurate qualitative and quantitative models— 
may be challenging to implement in unstructured envi¬ 
ronments or on imperfectly characterized machines. Sim¬ 
ilarly, methods depending on local linearizations of the 
typically (highly) nonlinear dynamics found in dynam¬ 
ically dexterous locomotion and manipulation systems 
[3, 4] typically suffer from small basins of attraction [5] 
and (to our knowledge) high sensitivity to parameters. ^ 

Observation (a) suggests that modularity of operation 
(i.e., wherein different combinations of actuators are used 
to effect distinctly different dynamical goals at different 
stages within the task cycle) will be a hallmark of practi¬ 
cal locomotion platforms. Observations (b) and (c) imply 
that simpler, less exact but potentially more robust repre¬ 
sentations of the principal dynamical effects likely to pre¬ 
vail across a wide range of substrates may offer a tractable 
means of working with rather than fighting against, or 

'in some robotics settings these disadvantages of the exact or local 
linearized control paradigm can be effectively remedied by recourse to 
parameter adaptation [6], but in our experience, such methods are too 
“laggy” to work in this hybrid dynamics domain with its intrinsically 
abrupt and rapidly switching characteristics. 






Figure 1: Control of a hopping behavior expressed as a hierarchical composition of closed-loop templates. Notionally, the grey 
arrows represent directed template^anchor relations. Center: A model of the tailed monoped physical platform on which we 
implement tail-energized planar hopping, labeled with configuration variables (black), actuators (red), and model parameters 
(blue). 


learning exactly the highly varied dynamical details. Ob¬ 
servation (d) implies that higher authority sensorimotor 
control activity ought to target continuous phases of the 
locomotion cycle, leaving the transition event interven¬ 
tions to more passive and mechanical sources of regula¬ 
tion [7]. In sum, these observations motivate the search 
for modular, reduced order representations of locomotion 
task constituents that are specialized to couple selected ac¬ 
tuation affordances to particular DOFs at particular phases 
of the locomotion cycle. The value of such component 
task representatives remains hostage to the availability of 
methods for composing them in a stable manner. 

This report introduces a novel locomotion platform, the 
Penn Jerboa, Fig. 7, to put a slowly maturing formalism 
for the composition of such modules to a practical test. 
We adopt the template-anchor^ framework [10] to repre¬ 
sent this machine’s 4DOF steady sagittal plane running as 
the hierarchical composition of the low DOF constituents 
depicted in Fig. 1 . At the leaves of this hierarchy tree, we 
introduce four different IDOF templates that emerge from 
the decades old bioinspired running literature [4, 11], 

^ The template-anchor relation as exemplied in various physi¬ 
cal [5, 8] and numerical [9] studies associates a pair of smooth vector 
fields, /^, on a pair of smooth spaces, “I c A via the condition 
that T is an attracting invariant submanifold of the anchor field, f^, 
whose restriction dynamics is conjugate to that of the template field, 

~ /'^ I a (where ~ denotes equivalence up to smooth change of co¬ 
ordinates). In this paper, we are dealing with hybrid fields and flows 
for which the extended definition and its verification is a bit more in¬ 
tricate. Thus exceeding the scope and length constraints of the present 
paper, we will treat the hybrid template-anchor relation as an intuitive 
notion here. 


joined by a new arrival from recent work on bioinspired 
tails [12, 13]. We apply the four decoupled IDOF con¬ 
trol laws associated with these isolated “leaf” templates 
directly to the (highly dynamically coupled) physical plat¬ 
form and demonstrate empirically steady sagittal plane 
running (on a circular boom) whose body motions reveal, 
when viewed in the appropriate coordinates, Fig. 15, strik¬ 
ing similarity to the corresponding isolated IDOF con¬ 
stituents. We show (up to a still unproven technical con¬ 
jecture) that the appropriate two pairs of these four IDOF 
leaf templates are formally anchored by the two “interior” 
2DOF templates depicted in Fig. 1, in the sense that the 
IDOF systems define attracting invariant submanifolds of 
the 2DOF systems that exhibit conjugate restriction dy¬ 
namics. We conjecture, as well, that the two interior nodes 
(the 2DOF templates) of the figure are in turn formally 
anchored by a physically realistic dynamical model of the 
closed loop Penn Jerboa in the sagittal plane. The data 
of Fig. 15 support this hypothesis, but we have not yet 
succeeded in completing the proof beyond the embedding 
and invariance properties. 

Notwithstanding the specifics of our compositional ap¬ 
proach to its control, we believe that the new physical 
platform is itself of independent interest by virtue of its 
added appendage (the “tail”), opening up a multiplicity of 
diverse uses for both of its two revolute actuators. Note 
again, however, this diversity of uses cannot be achieved 
without some recourse to behavioral modularity. In that 
light, we are particularly attracted by these simple low- 
DOF template controllers. In our experience, such con- 








Figure 2: Snapshots from apex to apex of tail-energized planar hopping (§5) implemented on a new robot platform—the Penn 
Jerboa (§6). 


structions have the hope of succeeding in unstructured 
outdoor settings, since they build on the relatively robust 
template dynamics. 

1.1 Relation to Prior Literature 

This “compositional” method of controller synthesis 
was pioneered empirically by Raibert [14] for planar and 
3D hopping machines, and we develop our planar hopping 
behavior by building up from those ideas. Our physical 
platform (Fig. 1 center) forgoes Raibert’s prismatic shank 
actuator, and instead places that actuator in an inertial ap¬ 
pendage. This motivates us to explore how tails can be 
“recycled” from their transitional agility duties [12, 13], 
now repurposed to substitute for Raibert’s shank actuator 
and play the role of steady-state running energizer in the 
sagittal plane. Apart from their use in transitional maneu¬ 
vers (inertial control in free-falling lizards [15] and robots 
[12,13] or in turning lizards [16] and robots [17]) it has re¬ 
cently been discovered that kanagaroos do positive work 
with their tails in a quasistatic pentapedal gait [18]. In 
our implementation, the tail contributes the reorientation 
function in flight, and the energetic “pump” function in 
stance (albeit in a dynamic fashion). We are not aware 
of prior robotic locomotion work wherein a tail is used to 
help power the stance phase. 

1.2 Contributions of the Paper 

This paper contributes both to the theory and practice 
of dynamical legged locomotion. 

The principal theoretical contributions are: (i) a new 
(slightly simplified) further abstraction (§3.3) of the long¬ 
standing SLIP running model [4] as a formal cross- 
product of previously proposed vertical [19] and fore-aft 
[20] templates; (ii) a stability proof (modulo a restrictive 
assumption 3) of the parallel composition^ of Raibert’s 
[14] stepping controller (10) with our new energy pump 
(3) in Proposition 6; and (iii) a proof of local stability in 

^By this term we mean the application to the (coupled) plant 
p^{x,u) (§3.3) of a decoupled control law, u = g'{xi) x g^^[x 2 ), 
taken directly from (3), (10), respectively. 


Table 1: List of Symbols 


i G 
D* 

/* : D* ^ TD* 
r* : dV* ^ 

F* : D* ^ aD* 

F* = f; o Fi* 
p*{x, u) 
h e 

a e 

R : ^ SO(2) 

Tx = {x, x) 

D^y 

fv G R+ 

G R ^ R + 

7 : R ^ 

/3 : R ^ 
ftw : ^ R^ 


Hybrid mode, where 1 is stance, 2 is flight 

Domain for template * in mode i 

Vector field in mode i 

Reset map from mode i to i -I- 1 

Mode i flow evaluated at the next transition 

Return map at touchdown (TD) event 

Plant to which we apply u = gi{x) to get /* 

Identity matrix of size d 

Planar skew-symmetric matrix 

T** standard basis vector 

Map from angle to rotation matrix 

Tangent vector associated with x 

Jacobian matrix dyijdxj 

SLIP radial velocity gain (§3.2.2) 

Map from radial TD velocity to k (§3.1.1) 
Fore-aft model stance sweep angle (§3.2.2) 
Raibert touchdown angle function (10) 
Cartesian to Polar TD velocity (§3.3.2) 


the inertial reorientation model (20) of the parallel com¬ 
position (21) of Raibert’s [14] pitch stabilizer and the tail 
reorientation controller [13] in Proposition 7. 

The empirical contributions of the paper are: (i) design 
and implementation of a working tailed biped platform, 
the Penn Jerboa (Fig. 7); (ii) physical demonstration of 
the (provably correct-Proposition 1) oscillatory spring- 
energization scheme for vertical hopping; and (iii) exper¬ 
imental evidence supporting the hypothesis that our final 
parallel composition of the four isolated controllers does 
indeed anchor the corresponding templates in the Jerboa 
body (Fig. 15). 

While the idea of parallel composition is appealing, 
the difficulty of such a composition arises from the nat¬ 
ural transfer of energy between different compartments 
[21]“^ in a mechanical system operating in a dynamical 

''We use this term here to stand for subsystems (here, disjoint sub- 
















Table 2: Template Controllers 


Tail energy pump 

g\{x) = kt cos(Zx) 

(3) 

Raiberl stepping [14] 

52^(x) = (3*{x) -F kp{x — i*) 

(10) 

Raibert pilch correction [14] 

~ kghCLX kg(l\ 

(21) 

Shape reorienlalion [13] 

5 ^( 02 , 02 ) = -kgka2 - kga2 

(21) 


regime. In our setting, some degree of coupling across 
compartments is crucial to the underlying design concept 
of driving the leg spring through torques generated “far 
away” in the tail. Thus, a naive approach of looking 
for exactly decoupled body dynamics is not fruitful^. In¬ 
stead, we analyze stability properties of (hybrid) closed- 
loop templates-which are not specifically associated to 
any body-without paying attention to the input structure. 
In agreement with intuition, we find (§5.4) fhaf minimiza- 
fion of cross-femplafe fransfer of energy-fhrough eifher 
fhe flows or fhe resef maps-resulfs in a successful compo¬ 
sition. 

2 Preliminaries: Organization and No¬ 
tation 

Table 1 confains a lisf of imporfanf symbols in fhis pa¬ 
per, including a sef of symbols for describing hybrid dy¬ 
namical sysfems. We adopf fhe modeling paradigm from 
Definifion 1 in [22], represenfing a hybrid dynamical sys- 
fem by fhe fuple (D, /, r) as defined in Table 1. We only 
consider fwo hybrid modes in fhis paper: ballistic flighf, 
and a sfance phase arising from a slicking conlacl af fhe 
“foe”. 

Superscripls on each of Ihese symbols denole fhe hy¬ 
brid template lhal if is a pari of, e.g. for confrolled 
vertical hopping (§3.1). The layouf of fhe paper roughly 
reflecls fhe lemplale-anchor hierarchy depicfed in Fig. 1. 
Namely, Ihere are fwo intermediate 2DOF lemplales—fhe 
SLIP, s, and fhe inertial reorienlalion, a—lhal comprise 
Ihe failed monoped, tm = {s, a}. They, in lum, are com¬ 
prised of fhe vertical, v, and fore-afl, fa, IDOF templates, 
s = {v, fa}, and respectively, Ihe shape, sh, and pilch, 
p, IDOF templates, a = {sh, p}. We endow Ihe IDOF 
templates al fhe lowesl level wilh an exemplar planl, wilh 

sets of the physical degrees of freedom) that exchange a resource (here, 
energy). 

^For instance, for hopping with the tailed monoped, the tail actu¬ 
ator and hip actuator seemingly work on differently “binned” tail and 
leg DOFs, but we energize the robot body with the tail through the leg 
spring. 


respecl lo which we will develop conlrollers for Ihe four 
template planls, in isolation. 

Sections 3-4 presenl Ihe 2DOF s, a templates lhal are 
direclly anchored in Ihe robol body (§5), and wilhin Ihem 
conlain descriptions of Ihe subtemplates (e.g. §3.1, 3.2) — 
as simple exemplar IDOF anchoring bodies and corre¬ 
sponding conlrol laws—lhal comprise in isolation Ihe 
consliluenl desired limiting behaviors lhal we seek lo em¬ 
body simullaneously in our physical system. Each of Ihe 
template conlrollers in Ihis suite is necessarily simple by 
dinl of ils origin as a feedback law for a highly abslracl 
IDOF lask exemplar. We hypolhesize lhal fhis combina¬ 
tion of algorilhmic simplicily and lask specialization may 
lend robuslness in Ihe empirical selling since conlrol poli¬ 
cies are nol sensitive lo, and certainly avoid cancellation 
of, forces arising from dynamical coupling in Ihe anchor¬ 
ing body. 

We emphasize lhal Ihese coupling-naive feedback laws 
(summarized in Table 2) are simply “played back” (mod¬ 
ulo scaling) in Ihe 6DOF body (§5) wilh all ils compli¬ 
cated Irue dynamical coupling. We show formally Ihrough 
various propositions in fhis paper lhal neverlheless Ihe sla- 
bilily of Ihe templates and subtemplates persisls Ihrough 
composition for Ihe dislal segmenls of Ihe free (Fig. 1) — 
SLIP as a composition of vertical hopping and fore-afl 
speed conlrol, and allilude slabilizalion as a composition 
of inertial reorienlalion and Raibert’s pilch conlrol. We 
provide some preliminary suggestions aboul Ihe composi¬ 
tion of SLIP (s) wilh allilude (a) comparlmenls (center of 
Pig. 1), bul a full analysis is left lo fulure work. However, 
we offer empirical dala in §6 showing how Ihis idea has 
resulted in promising qualilalive behavior on Ihe Jerboa 
robol (Pig. 15, video allachmenl). 

3 The (2DOF) SLIP Template 

3.1 Controlled Vertical Hopping (IDOF) 

Lor a successful hopping behavior, energy musl be 
periodically injected into Ihe robol body to compensate 
for losses. We simplify Ihe analysis here to a IDOL 
verlically-conslrained poinl-mass which can alternate be- 





Table 3: Physical Parameters (all scalars unless noted) 


kt Tail gain (3) 

kp Raibert speed controller gain (10) 

k Inertial reorientation generalized damper gains (21) 

kg Inertial reorientation graph error gain (21) 

(T, oj Dissipation, frequency of spring-damper (§3.1) 

£ Saturation parameter for tail controller (3)) 

Sr Stability margin for vertical hopping (Proposition 6) 
Sa Arbitrarily small orientation error (Proposition 7) 
mb, ib Mass, inertia of robot body (§5) 
pi,pt Leg, tail link lengths (§3,5) 
ks Hooke’s law leg spring constant (§3,5) 


tween stance phase (during which the actuator has affor- 
dance) and a ballistic (passive) flight phase. It has been 
shown in the past empirically [14] and analytically [23] 
that an impulse at the bottom of stance can produce a sta¬ 
ble limit cycle, in the presence of a spring for energy stor¬ 
age. In this paper, we consider a different strategy of an 
actuator forcing the damped spring by applying forces in a 
phase-locked manner. This choice of input representative 
is made with an eye toward using a tail actuator exerting 
inertial reaction forces on the spring (this model is for¬ 
mally instantiated §5). Intuitively, this can be thought of 
as negative damping [19] (effectively cancelling losses by 
physical damping). 

Throughout this paper, we make the following assump¬ 
tion inspired by [14]: 

Assumption 1 (Stance duration). The duration of stance, 
Tg, is approximately constant. 

This essentially asserts that the damping losses or ac¬ 
tuator forces are relatively small compared to the spring- 
mass dynamics (in their effect on the liftoff condition). 

We build upon the “linear spring” analysis in [23] for 
our vertical hopping exemplar body and closed-loop tem¬ 
plate. For a spring-mass-damper system with spring de¬ 
flection Xj damping coefficient ^ and natural frequency 
u 

X + 2w/3x + = r. (1) 

With the change of coordinates xi := X’ ^2 '■= xju, 

X = r) := —ujJx + {— 2 ^ 00 x 2 + Tjuf), (2) 

and the hybrid reset events occur at xi = 0 (corre¬ 
sponding physically to the touchdown and liftoff events 

at X = 0). 



Figure 3: Left: The vector field and an execution of (4), show¬ 
ing a stable limit cycle. Right: The vertical “energy” is easy to 
tune with kt ■ 


3.1.1 Oscillatory Spring Energization 

We choose the physically motivated control strategy 

T ■= ^ hcos/^x, (3) 

where e > 0 is a small saturation constant. It is clear 
in this form that the input is a fed-back version of the 
“phase” only. We obtain the closed-loop stance dynam¬ 
ics 


i = fi (x) ■■= -ooJx + (^-2/3uj + X2e2. (4) 

Proposition 1 (Oscillatory energization stability). The 
vertical hopping template (4) has a unique attracting pe¬ 
riodic orbit. 

Proof. First, note that a: = 0 is the only equilibrium 
of (4). Secondly, note that 

x^i = xl (-2Pu -f , (5) 

which is zero on the set ||xj|* = — e. Additionally, 

since x'^i:||| 3 ,||<||a,||* > 0 and a:^i|||j,||>||a,||* < 0, this limit 
cycle is attracting. □ 

Writing x{t, xq) to denote the flow generated by (1), 
and letting S(xo) := min{f > 0 | eJx{t,xo) = 0} denote 
the stance time (since x\, vanishes exactly at the liftoff), 
we define the vertical stance map, 

-^1 (x) := ei’x(S(x, 0), (x, 0)). (6) 

As a corollary to Proposition 1, we know has 
an asymptotically stable fixed point, x* , and —1 < 
< 1 . 

Ballistic flight simply reverses the velocity, 

Ff{x) := -X- 


(V) 













Note that by symmetry (f^, and consequently are 
odd), o F{ = F^ o F{ o F^ o F^, i.e. the stability 
properties of the hybrid system are the same as that of the 
stance map as analyzed in Proposition 1 . Define 

« = Kix) ■= ( 8 ) 

the effective coefficient of restitution through stance, or 
the so-called “velocity gain” during SLIP stance [20]. 
Note that there is a unique fixed point, k* = 1, in these co¬ 
ordinates, which is necessary and sufficient for the smooth 
invertibility of h^, as can be seen by direct computation 
of its derivative. 

Conjugating the touchdown velocity return map via this 
diffeomorphism, we can define a return map for n, F'', 

F''{k) := o F^ o F{ o h~^{K) = h^{Kh~^{K)). (9) 

Proposition 2 (Vertical stability). The velocity gain re¬ 
turn map, F'', has an asymptotically stable fixed point, 
K* := 1, and DF''\f^^i = —DF^\j^*. 

Proof. This directly follows from the observation that n 
and touchdown velocity are related by a diffeo. Proposi¬ 
tion 1, and the simple form of F 2 in (1). □ 

3.2 Controlled Fore-Aft Speed (IDOF) 

Running and walking systems of a large variety from 
the sagittal or frontal plane resemble inverted pendula dur¬ 
ing stance [4], usually controlled by stepping strategies. It 
has been shown that a fixed touchdown angle can admit a 
reasonable basin of stability around an emergent attract¬ 
ing steady-state velocity in SLIP [24] . The capture point 
[25] and zero moment point [26] methods use a quasistatic 
heuristic which is related to these ideas, but are not ex¬ 
plicitly designed to servo to desired nonzero speeds. We 
attempt here to place the empirical success of [14] in the 
context of a model where its stability properties can be 
analyzed. 

3.2.1 The Raibert Stepping Controller 

In his classical empirical study, Raibert [14] inspired 
decades of subsequent experimentation and analysis by 
offering the following observations® about the pendular 
stance phase in his running machine travelling at forward 
speed, X, and stepping with a touchdown angle j5{x) (as 
in Fig. 4): 

Assumption 2 (Raibert observations). (i) For each speed, 
X, there is a neutraf touchdown angle, fi*{x) (ii) this 

^These conditions are not a direct result of SLIP’s nonlinear dy¬ 
namics, but are applicable to regime of interest. 

’in this context, “neutral” means = x, where x'^ refers to the 
fore-aft speed at the subsequent touchdown event. 



Figure 4; A simple model for the IDOF fore-aft dynamics in 
SLIP, closely related to BHop [20]. 

neutral angle is monotonic with speed, > 0, and 

(Hi) deviations from touchdown angle cause negative ac¬ 
celeration, i.e. — x)\j 3 ^i 3 * < 0. 

Proposition 3 (Raibert stepping controller). Under as¬ 
sumptions 2(i-iii), the Raibert stepping controller, 

/3 : X >->■ I3*{x) kp{x — x*) (10) 

stabilizes the forward speed to x*. 

Proof. Note that 

D±(x'^ - x) = Dp{x'^ - x) ■ D±fi{x) 

= - x) ■ + kp) 

= 1 -f - x) ■ (Difi* -t- kp). 

From the sign properties of various terms, we note that for 
small kp, —1 < D±x^ <1. □ 

3.2.2 Modified BHop as a Fore-Aft Model 

Building on existing SLIP literature [27], we make the 
following assumptions about pendular stance: 

Assumption 3 (Pendular stance). During stance, (i) the 
effects of gravity are negligible^ compared to spring po¬ 
tential / damping forces, (ii) radial deflections are negli¬ 
gible, (Hi) time of stance is constant, and (iv) the angle 
swept by the leg admits a small-angle approximation. 

Schwind [27] approximated that angular momentum 
about the toe is constant during stance, but we simplify 
further with the second assumption, and conclude that the 
angular velocity is roughly constant during stance. We 
adopt the third approximation from Raibert [14], and the 
last approximation is made for the ensuing analytical sim¬ 
plifications in §5.4, but we find empirically (§6) that it is 
not critical in practice. 

*We suspect that the less restrictive Geyer approximation [28] is 
sufficient, but leave this generalization to future work. 










These assumptions lead directly to the construction of 
the following return map acting on touchdown velocity in 
Cartesian coordinates (cf. Fig. 4). Then, 

F%v,k) = [^_,] R(-7 + /3) [ 1 _.] R(-/3)u 

= ( 11 ) 

where k (explicitly, the interaction from the radial compo¬ 
nent of SLIP) is taken to be a fixed parameter at this stage, 
7 (^ 1 ) is angle swept by the leg over the course 

of stance and /3{vi) is the leg touchdown angle (§3.2.1). 
This model is only a slight modification^ of BHop [20]. 

This analytically tractable model (i) allows us to “sep¬ 
arate” the radial dynamics (encapsulated in k) from the 
contributions of the fore-aft model itself, (ii) captures the 
exchange of vertical and horizontal energy through step¬ 
ping, and (iii) matches the empirically observed Raibert 
conditions (Fig. 5) as well as empirical data (Fig. 15), 
suggesting it is physically applicable and not just an ana¬ 
lytical convenience. 

For now we restrict our attention to «: = 1, and gen¬ 
eralize to include the radial dynamics in §3. With this 
restriction. 



Figure 5: A contour plot of the fore-aft acceleration a;+ — i 
produced by the MBHop model for a range of fore-aft speed 
X and touchdown angle (3. This plot depicts that (in a range 
around the neutral angle), this model captures all the conditions 
of assumption 2. 


which is lower-triangular. The eigenvalues are {1,1-1- 
2/cpu|}, which capture the local stability of the single 
fore-aft DOF (1 -I- 2kpV2 < 1) as well as the degeneracy 
of the map. 

To see why the last statement is true, note that we can 
find a rank 1 map 


:=F^(u,l) =R(7-2/3 )u, (12) 


i ^ M : u i-> !|u 


While we choose fo parameterize fhe refurn map as a 
funclion of u 6 M^, if is really a ID map: 


which is invarianf fo i.e. l = loF^^. Taking a gradienf 
of bofh sides and using fhe chain rule. 


Proposition 4 (Fore-afl sfabilify). MBHop with the Raib¬ 
ert controller presents a stable touchdown return map. 

Proof. We can check fhaf F^^ safisfies each of fhe Raib¬ 
ert conditions (Fig. 5), fhereby concluding aufomafically 
from Proposifion 3 fhaf fhe Raiberf confroller will ensure 
local sfabilify. 

Alfemafively, fhe ufilify of our simple analytical model 
( 11 )-( 12 ) is fhaf we can direcfly compute fhe sfabilify 
properfies under fhe Raiberf confroller (10), 

DF^^{v) := R + JRv • (D 7 - 2D/3)ef, (13) 


where R is evaluafed af 7 — 2/3. By inspection, fhe (de¬ 
sired) fixed poinf of (12) is /3 = 7/2 (fhis is fhe neufral 
touchdown angle). Evaluafed af fhe fixed poinf. 


DF^^v*) = I-2kpJv*el 


l+2kpV^ 0 " 
- 2 kpV* ij ’ 


(14) 


® Specifically, the similarities are apparent between (11) and (19) of 
[20]. The slightly discrepancy should be attributed to our insistence on 
using the physical touchdown and sweep angles /3 and 7 in the model, 
whereas the abstract parameter 6 in [ 20 ] results in a more succinct 
form. 


Dily — ■ DF^^\y. 

Evaluating at the fixed poinf v*, 

= F)t^y^ ‘ F)F 

i.e. Di\y* is a leff eigenvector of DF^^\y* wifh unify 
eigenvalue. 

Consequenfly, under iferafions of fhis map, we gef an 
invarianf submanifold spanned by fhe orthogonal comple- 
menf of fhe unify eigenvecfor, resulting in a “dimension 
reducfion” (fo a codimension 1 submanifold). In our case, 
F^^ is really a ID map, even fhough ifs (co)domain in M^. 

□ 


3.3 SLIP as a Parallel Composition 

In order to anchor our IDOE templates in the classical 
SEIP model (2DOP point mass with 2DOP springy leg), 
we simply “play back” our devised control schemes (Sec¬ 
tions 3.1 and 3.2). In the following subsections, we check 
that the closed-loop executions in the higher-DOE body 
still resemble a cross-product of our template behaviors. 





















For instance, prior literature has observed a decomposi¬ 
tion of SLIP dynamics into radial and tangential compo¬ 
nents, but to our knowledge there is no complete account 
of the stability of the parallelly composed (closed-loop) 
templates in these components. 


3.3.1 Hybrid Dynamical Model of SLIP 

We will construct our template plant model from [27] : 
a bead of mass 1 at (Cartesian) coordinates {x^, z^) e M^, 
with a springy (Hooke’s law spring constant kg) massless 
leg of length*® e M+ (where M+ is restricted to strictly 
positive reals, and is open) and rest length pi, at an angle 
of e S'* from vertical. Let := (0f, x^, z^). Using 

assumption 3(iv) as a convenience (though that assump¬ 
tion is not required for this formulation), the touchdown 
and lift-off conditions can be specified in terms of the ze¬ 
ros of a® := z^ — pi- 

Define Q| := S* x M_|_ x M x Jj, where M = Ji u J 2 := 
{-CO, pi] u {pi, co). Then, D- := TQ-, and 


/f(q=,q^) 



eiq +ks{pi-ei) 


(15) 


/Kq'iq") 



(16) 


where fhe unspecified componenfs are (i) fhe mass-cenfer 
dynamics which are consfrained by [ ] = 01 

in (15), and (ii) fhe degenerate massless leg dynamics in 
(16). 


The Guard Set is dD^ Since is ilself a cross producf 
of Euclidean spaces and Lie groups, we can identify fhe 
fangenf bundle wifh a cross producf, TQ^ Q| x M^. 
Then, fhe boundary of fhe product space only contains 
parts from Jj, which corresponds exactly to the zeros of 
aM§3.3.1). 


Reset Maps Let us define fhe functions 


Cart ; 5* x 


Pol ; 


02 


— sin 9i 
cos 9i 


5* X 


U I 


Zit 

ll“ll 


(17) 

(18) 


'®We use 0 for leg “joints” to be consistent with [29]. 


The resef maps are defined as 


: D" 




e 

e 

[S] 

[i] 


D" ^ 


e 

e 

["] 

[i] 


9 

9 

Cart(0) ; 

_ DCaTt\0-9 _ 

DPol-[|j 
1^ — 2 : tan /3(i) j 

. [ 5 ] 


3.3.2 Anchoring the IDOF Templates 

Consequent upon the above model—where each hy¬ 
brid mode is dynamically 2DOF—SLIP is a 4D dynam¬ 
ical system (one parameterization being {x,z,v), where 
n 6 is the touchdown velocity, and {x, z) e is 
the Cartesian location of the point mass at touchdown). 
The efficacy of our 2D refurn map analysis is esfablished 
by argumenfs similar fo fhose of [30]: fhe Poincare sec- 
lion z^^ = Pi cos j5{v) eliminates one dimension, and fhe 
equivariance of fhe dynamics wifh x eliminafes anofher. 

We firsl observe fhaf our MBHop model of §3.2.2 still 
represents the pendular stance correctly under assumption 
3. However, k is not a fixed parameter, bul evolves ac¬ 
cording fo dynamics similar fo F'^ in Proposifion 2. From 
(10) and (11), fhe embedded («: = !,?; = r;*) submani¬ 
fold is invarianf. We show in Proposifion 6 lhal if is also 
aflracling. 

Fef us define ^ as 


w = hvj{v) := R{—/3{v))v. (19) 

Lemma 5. Let V := {n 6 : r ;2 < —Then hw\v R 

a local diffeomorphism. * * 

Proof. Nofe lhal 

Dhvj = R — JRvDPci, 

where R is understood to be evaluated al —/3{v). By in¬ 
spection, 79hw could only have a lesl vector R^ JRn in ifs 
kernel, i.e. 

Dh^ ■ (R^JRr;) = (1 — D/3eiR^JRv)JRv 0, 
since we know v ^ 0, Djd = and so 

1 — 79/3efR^JRr; = 1 + V 2 + kp^ < 0, 

by fhe conditions assumed on kp. Thus 79/i^ is nonsingu¬ 
lar, and hw is a local diffeo. □ 

"Physically, the restriction to V means that the hopper must have 
sufficient vertical component of touchdown velocity, essentially elim¬ 
inating “grazing” ground impacts. 

























The vector w gives a tangential/radial decomposition of 
V (i.e. polar with respect to the leg angle). 

Additionally, using (8), we can “recover” the k- 
dynamics in the coupled system: k = hi^{w 2 )- We prefer 
the redundant {v, k) parameterization because of analyti¬ 
cal tractability. 


Now notice that since Dh^ is well-conditioned, we can 
claim an upper bound on 

\p^q\ ^ 2/cpll ^ kpE. 


Also, the quadratic form 




p must have 


Proposition 6 (Stability of SLIP as a composition). For 
(i) stable vertical hopping with — 1 -I- er < — < 

1 — ffr. (ii) sufficiently^^ small kp in the Raibert contoller, 
parallel composition of the radial and fore-aft templates 
results in a locally stable 2D return map, F^. 

Proof We choose to perform our stability analysis at a 
section just after touchdown (in w = h,f,{v) coordinates). 
From (11), the return map in m-coordinates is 


T 

" f 


q 

-DF-\-\ 

P 




since 


has norm less than 1. 


It can be checked that both eigenvalues are of absolute 
value bounded by unity iff all of (i) det < 1, (ii) det > 
tr — 1, and (iii) det > —tr — 1 are true. These inequalities 
follow from condition (ii) of Proposition 6 and choosing 


small enough kp such that 2kpE < 8^ 


□ 


F^{w) := /iw o o /iw iw)\K=h 4 w 2 ) 




hK.{w2) 


W, 


where t] ■.= [y — (5 — j3 o F^) o Now, 


DF^ = D^F^ + D^F^ • 


where the first summand can be thought of as loosely the 
isolated fore-aft subsystem behavior, and the second sum¬ 
mand is the perturbation from the radial subsystem. We 
will evaluate this quantity at the fixed point = hy„{v*). 

Observe that using (10), Drj\^ = —2kpe{Dh~^. Pro¬ 
ceeding just like in Proposition 2, 

Dhf\* = L|) > 

Df^F^ = R(?7)e2ejt(; D^F% = wle2. 


Lastly, the “isolated” term computes similar to (14), 

D^F^ = R[i ^] + JR[i 

DwF^\^ = I + Jw*Dp\^. 


Putting all of these together. 






+ pq^, 


where p := —2kpJw*, q^ := ejDh,^^. Using the matrix 
determinant lemma. 


tr DF^ = 1 — DFf I * + p'^q 


det DF^ = -DFf\,, (l 


1 

-DF-\-^ 



'^Formally, this means that kp can be chosen as a function of Er- 


4 Hybrid Inertial Reorientation 
(2DOF) 

Our decision to energize the hopping behavior with a 
tail leaves introduces a new actuated DOF whose tight 
dynamical coupling to both the mass center and the body 
orientation dynamics requires its careful control through¬ 
out the locomotion cycle. Recent literature [13] has seen 
the development of a IDOF “inertial reorientation” tem¬ 
plate for correcting the “shape” coordinate in a two-link 
body experiencing free-fall (constrained by conservation 
of angular momentum). Raibert [14] introduced a pitch 
stabilization mechanism relying on reaction torques from 
hip actuation during stance. In this paper, we adopt the 
approach of composing these templates for 2DOF stabi¬ 
lization of appropriately defined “pitch” and “shape” co¬ 
ordinates of a two-link body/tail model. 

Since in the physical system the tail actuator, T 2 , is un¬ 
available for attitude control in stance (because it is being 
“monopolized” as the destabilizing energy source for the 
SLIP subsystem), and the Raibert pitch correction mech¬ 
anism (using the hip actuator, rf) is unavailable in flight 
(due to absence of ground reaction force), we present a 
hybrid inertial reorientation (HIR) template (Fig. 6) as the 
simplest exemplar body on which this 2DOF template is 
anchored. 

We omit the Lagrangian derivation for this familiar sub¬ 
system [13], but exploit the fact that when pinned at the 
CoM, the dynamics are second-order LTI with no Coriolis 
terms. We perform a change of coordinates (inverting the 
constant inertia tensor) to obtain the (decoupled) dynam¬ 
ics 




ai 

= J 

a2 

) 




[7] =-P\{Ta,Ti) (stance), 

[Ti] =-Pl{Ta,T2) (flight), 


( 20 ) 



















Figure 6: A hybrid 2DOF inertial reorientation template with two segments pinned at the CoM and no gravity. Left: the net 
angular momentum of the system is constant. Right: the system can correct the net angular momentum using reaction torques 
on the main body segment, but the tail DOF is subject to an unmodeled disturbance , or (5 in (20). 


where (01,02) are the “pitch” and “shape” coordinates, 
respectively, and 6 is an unmodeled disturbance term (ex¬ 
plicitly added here with an eye toward the use of tail for 
spring energization in the physical system). In (20) we 
have now represented HIR as two independent subsys¬ 
tems on which two identical IDOF templates will be an¬ 
chored in parallel (albeit in alternating stages of the hybrid 
execution). 

Taking advantage of the direct affordance (by which 
we mean that both of the two decoupled IDOF systems 
are completely actuated in, one and then other, of the al¬ 
ternating modes of their hybrid dynamics), we employ a 
graph-error controller [31] as a type of reduction. Since 
our reference first-order dynamics are just di = —kai, the 
independent closed-loop IDOF subtemplate vector fields, 
/P : Tai t—> Toi and : Ta 2 Ta 2 , are defined as 


di = —kg{di + ktti) = —kgktti — kgdi, ( 21 ) 


where the gain kg is understood to be high enough to make 
the transients of the anchoring dynamics irrelevant. 

4.1 Hybrid Dynamical Model of HIR 

Since the isolated model does not have any intrin¬ 
sic physical mechanism for transitioning between modes, 
we add an exogenous clock signal, V’a £ such that 
'i/’a £ [0, vr] represents stance, and the complement repre¬ 
sents flight. In this paper we sidestep the issue of phase- 
synchronization for the various compartments, but sim¬ 
ply use ipa to ensure our gains our tuned properly for the 
timescales of the coupled system (Proposition 9). 

Define = TS^x {(0, vr] u (vr, 27r]}. Now the closed- 
loop template dynamics, : TS'^ x ^ T{TS^ x S^) 


can be specified as 
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( 22 ) 


the guards sets are dV^ = TS'^ x {{tt} u {27r}} and the 
reset maps r| = id simply modify the dynamics (20) at 
ij^a = (stance to flight) and '0a = 0 (flight to stance). 


4.2 HIR Stability Analysis 

Let us denote := J ddt, the interval being over the 
stance phase of stride i. Also, define 5max = maxt 5[f]. 

Proposition 7 (HIR Stability). Setting 

> 2^ log (1 + 6 m^^/ea) 

results in the desired limiting behavior for F^: ||o|| —> 
(0), a neighborhood ofO of size Sa- 

Proof Simply integrating the first-order dynamics (22), 
we get the touchdown return map : 3“^ ^ 

F%a) = C-{a + 5[^,]), (23) 

where ( := — kiTjuia)- Iterating this return 

map, at stride n 6 Z+, 

a[n] = C*a[0] -f (^^[O] H--f C6[n - l])e2, (24) 


and using the triangle inequality. 


«MII ^ ICr • l|o[0]|| -f (5max 


(25) 


Note that C < , , ? ^ ^— is a sufficient condition to en- 

l + Dmax/^a 

sure that ||a[f]|| ^ Sa asymptotically stable. Some algebra 
reveals that 


k > 



(26) 


is, in turn, a condition sufficient to insure that previous 
inequality involving □ 
















Tailed monoped 

Figure 7: The Jerboa is a 2 Kg robot with hip-actuated legs 
and a 2DOF tail, pictured on the left as it appeared in the exper¬ 
iments of Section 6. On the right is our model for the planarized 
4DOF system for comparison. 


5 Physical System: Tailed, Compliant¬ 
legged Biped 


Our target physical platform is a tailed bipedal robot 
that we have built, which (when planarized) we model 
as shown in the center of Fig. 1. We were able to for¬ 
mally show template-anchor relations going from IDOF 
to 2DOF templates (Propositions 6 and 7), because of the 
availability of simple models (§3.2.2), or trivial dynamics 
(§4). However, as we proceed up the desired hierarchy 
(Fig. 1), there are no easily accessible tools that let us di¬ 
rectly analyze the effects of coupling in the return map. 
In this section, we only show (Proposition 9) that under 
a highly restrictive assumption 5 (that essentially makes 
the tail sweep negligible), the closed-loop tailed monoped 
return map has an invariant submanifold where it is 
equal to x F^, but we also leave as conjecture that this 
invariant submanifold is attracting. 

The first two subsections of this section discuss (in an 
informal manner) the design process of the robot plat¬ 
form we have designed, built, and implemented the tail- 
energized hopping behavior on. 


5.1 Jerboa: Design and Construction 

The Jerboa was designed with the goal of being a dy¬ 
namic, agile robot with an inertial appendage. We defer 
an in-depth discussion of morphological constraints and 
tradeoffs to future work, but present the following basic 
design decisions here: 

i) With an eye on power density constraints'^, the robot 
is underactuated. There are 12 spatial DOFs (6 for 
the body, 2 for each re volute leg, 2 for the tail) and 
4 actuators. When planarized with a boom and vir¬ 
tual constrains on the appendages (as we have done in 
this paper), there are 6 planar DOFs: 3 for the body, 
2 for the single leg and 1 for the pitching tail. Raibert 
showed that an underactuated robot can be dynami¬ 
cally stable [14], and in order to have the best per¬ 
formance, we limited the number of actuators on the 
robot to the minimum number that we believe is re¬ 
quired to achieve a wide variety'^ of behaviors. 

ii) The body has low inertia (due to the mass of the 
motors being concentrated near the CoM, and the 
appendages being light), and the actuators are con¬ 
figured such that they can impart correspondingly 
large accelerations to the body (with an eye to¬ 
wards “agility”). Future work is planned to reconcile 
our inclination with emerging definitions of specific 
agility [32], but intuitively it seems as if “integrated 
magnitude of body acceleration” is a reasonable met¬ 
ric to aim for. 

iii) The hips are actuated, but the leg extension is com¬ 
pletely passive. This particular form of underactu¬ 
ated leg has been demonstrated to have great versa¬ 
tility in RHex [33], for steady-state running as well 
as transitional maneuvers [34] . 

iv) The robot contains an inertial appendage which is 
endowed with the same amount of power as the 
hips. Recent research in biomechanics [12] and 
robotics [13] has demonstrated the utility of tails as 
inertial “self-righting” devices, and on the Jerboa we 
promote it to a primary source of locomotory energy 
and control. 

In the remainder of this section, we outline the elec¬ 
tromechanical aspects of the construction of the robot. A 

'^Adding actuated DOFs parasitically increases mass, but it is not a 
direct consequence that a proportionate amount of usable power will 
be added to the robot body by the extra DOFs. 

'"'We have some preliminary empirical evidence that the Jerboa can 
quasistatically and dynamically balance, in order to sit, stand, walk, 
hop, run, turn, leap, etc. Careful investigation of each of these behav¬ 
iors is planned for future work. 






Figure 8: The jerboa tail is a 2DOF spherical joint controlled 
using coaxial motors through a mechanical linkage. Though 
there are driven sprockets visible in this image, the version of 
the robot presented in this paper did not have this additional 
reduction stage. 

summary of important mechanical measurements is pro¬ 
vided in table 4. 

5.1.1 2DOF tail 

The tail appendage is configured as a 2DOF spherical 
joint with a point mass at the distal tip. The joint itself 
is constructed using a linkage (Fig. 8) such that identical 
motor displacements result in a pitching motion, and dif¬ 
ferential motor displacements result in a yawing motion. 
The forward kinematics map from motor angles fti, /i 2 £ 
to the tail pitch and yaw angles, ( 1 ) 2 , 0yaw e has a 
simple form when restricted to zero yaw (i.e. /ii = ^ 2 ), 

Ai2)Uya,w=o = fti = F2- (27) 

For the behavior under study in this paper, a virtual con¬ 
straint ensures that (/>yaw = 0. We leave a full kinematic 
analysis of the 2DOF mechanism to future work. 

5.1.2 Prismatic-compliant revolute-actuated legs 

Even though we adopt the underactuated hip-driven 
legs from RHex, the legs are chosen to have “toes” with 
point contacts instead of the rolling contact typical of 
RHex legs for the following reasons: 

i) Our template plant for fore-aft speed control (§3.2.2) 
is an inverted pendulum with a point contact, and 
in particular, the toe-placement strategy for fore- 
aft speed control (10) is only (currently) well- 
understood for this leg structure. 

ii) The Raibert pitch controller [14], which we use as 
part of our attitude control (21), depends on a “rigid” 
connection between the hip and the toe. With a 
series-elastic element that may have torsional com¬ 
pliance (such as a C-leg), the ground reaction force 
would load up the leg spring, introducing the spring 
dynamics as a “lag” in our pitch control strategy. 


The left of Fig. 9 shows three leg designs that were con¬ 
sidered for the Jerboa: i) a prismatic mechanism with a 
nonlinear elastic element, ii) a compression spring in a 
four-bar mechanism, and iii) an extension spring in a four- 
bar mechanism. While the kinematic properties of the first 
design are the closest to our model (the spring force at 
the toe is purely radial, and the motor force at the toe is 
purely tangential), this design proved difficult to construct 
because of the linear bearing required. The kinematics of 
the “approximate” leg designs are pictorially depicted in 
Fig. 9. 

The experiments for this paper were all performed 
with the compression-spring legs. The compliant element 
is an off-the-shelf shock absorber for RC vehicles with 
lightweight construction, but considerable damping. We 
believe that the damping in legs was an important limiting 
factor in the energy of the hopping behavior demonstrated 
in §6. 

5.1.3 Actuators 

The power generated by electromechanical actuators 
tends to be at unusably high speeds for legged applica¬ 
tions, however at the same time, higher gear reductions 
are undesirable due to a multitude of reasons [35]. To this 
end, we tune our actuator selection to maximize thermal 
specific torque, TFts— the torque generated by the motor 
per unit mass per unit temperature rise. This modification 
to the torque density criterion of [35] allows us to incor¬ 
porate the thermal implications of sustained motor activa- 
tion^^. Fig. 10 contains a table comparing these metrics 
for the chosen actuator, a T-motor U8, and the one se¬ 
lected for X-RHex [36], a Maxon EC-45. 

Additionally, we developed custom motor con¬ 
trollers built around Infineon BTN8980 integrated half¬ 
bridges and an STM32E373 microcontroller that are 
(a) lightweight (20 g), (b) commutate using field-oriented 
control (EOC) at 25 KHz (adapted from [37]), (c) deliver 
up to 55 A peak current and up to 40 V peak voltage, and 
(d) have built-in 12-bit rotor position sensing. As a trade¬ 
off for the high power-density of the driving electronics, 
they are limited by the heat dissipation ability of the half¬ 
bridges. Based on some crude testing, we have found that 
we can source approximately 10 A of steady-state cur¬ 
rent (thermally limited), corresponding to around 1 N-m 
of torque. Eig. 11 compares the physical dimensions and 
thermal performance of the motor controllers to the mo¬ 
tors we have chosen. We note the following consequences 
of our selection of motor and driving electronics: 

'^We are assuming a thermal dissipation model for the motor, but 
not accounting for temperature effects on magnetic flux density. 






Table 4: Parameter values 


Mass (with battery) 

2.419 Kg 

Dimensions (without tail) 

0.21 m (L) X 0.23 m (W) x 0.1 m (H) 

Tail length 

0.3 m 

Tail mass 

150 g 

Leg length 

0.105 m 

Leg motor stall torque 

3.5 N-m 

Peak power density 

376 W/Kg 

Peak (vertical) force density 

46 N/Kg 



Figure 9: Left: Three leg designs considered for the Jerboa; the prismatic spring is “ideal” (in our model of §3, the spring force 
is dominantly axial, and the actuator force is predominantly tangential) but difficult to manufacture, and the four-bar designs 
only approximate the desired kinematics. Right: Configuration-dependent Jacobians of the compression and extension spring 
designs, where the displayed arrows map infinitesimal hip torques and spring extension forces to forces represented by red and 
blue (resp.) arrows at the toe. Out of these designs, the pictured version of the robot in Fig. 7 uses compression springs. 


i) The high torque density of the chosen motors allows 
us to completely forgo any static gear-reduction on 
the Jerboa (although the 2DOF tail makes use of 
a linkage to transmit power to a spherical joint)— 
affording benefits of “transparency” and eliminating 
any transmission losses [35, 38]. 

ii) Power dissipation (to heat) in the motors is not a lim¬ 
iting factor in the robot’s performance with the cur¬ 
rent driving electronics. 

hi) By eliminating the need for gearboxes and judicious 
chassis design, we have been able to reduce the 
“robot framing cost” to only 40% of the mass of the 
robot. To put this in context, only 8% of the mass of 
X-RHex is motors [36]. 

Lastly, we highlight some of the design aspects of the 
Jerboa that are particularly relevant to the subject of this 
report (tail-energized hopping via parallel playback of de¬ 
coupled controllers): 

Assumption 4 (Design for decoupled control). The de¬ 
sign of the Jerboa specifically ensures (i) leg/tail axes 


of rotation are coincident at the “hip,” (ii) the tail mass 
is small, i.e. mt « nib, and (Hi) center of mass (ap¬ 
proximately configuration-independent by the previous 
assumption) coincides with the hip. 

We point out here that these design decisions are less 
strict than the ones required for our present analysis (as¬ 
sumption 5). We believe that the stringency of assump¬ 
tion 5 is not necessary, and provide some empirical evi¬ 
dence to this effect in Section 6. 

5.2 Modeling for Planar Hopping 

Raibert’s planar hopper [14] empirically demonstrated 
stable hopping using a rigid body with a springy leg, and 
in this paper we pursue the same idea, but instantiate verti¬ 
cal hopping by coupling the 1-DOF leg-spring excitation 
controller (physically acting through the tail). In flight, 
the tail actuator grants us a new affordance that we only 
use here to regulate the added “shape” DOF. Our phys¬ 
ical model is shown in Fig. 7. The system has a sin¬ 
gle massless leg with joints 0 = { 61 , 62 ) e x M+, 
a rigid body {x,z,(j)i) e SE(2), and a point-mass tail 
with revolute DOF (/) 2 , such that the full configuration is 







































Figure 10: Left: The selected actuator for the Jerboa is the T-motor U8, showing a thin profile and large gap radius—desireable 
properties for legged applications [35]. Middle: Motor properties relevant to selection for legged applications for the Jerboa 
motor, and the X-RHex [36] motor. Right: A torque-current plot for the U8 when coupled with our custom motor controllers of 
Fig. 11, showing flux saturation at higher currents and a dashed line for the nominal torque (predicted by Kt)- 



Figure 11: Left: The physical dimensions of our motor controller when compared to the motors they are driving. Right: Infrared 
image of of our actuation setup at stall, showing the controller reaching higher temperatures than the motor coils. 


q := { 6 i, 92 ,x,z,(j)i,(j) 2 ) e Q- We make the following 
design-time assumptions: 

Assumption 5. (i) Leg/tail axes of rotation are coincident 
at the “hip,” (ii) tail mass is small, i.e. nit « w-fo. (Hi) cen¬ 
ter of mass (configuration-independent by the previous as¬ 
sumption) coincides with the hip, and (iv) body, tail have 
high inertia, i.e. ib,it —> oo. 


such that Ai(q) = [RDg/JRgo], In flight mode, 
a 2 {q) = 0. As in [29], the dynamics can be expressed 
as 


^ 1 [q] =[T-N] 

A, 0 J LaJ L 0 J 



(30) 


Define the linear coordinate change h : y = S x A —> 
Q, and H := Dh such that 


5.3 Equations of Motion 

Using the self-manipulation [29] formulation of hybrid 
dynamics, the inertia tensor is 


M = [°mJ > where Mb 


Ml 

Mo M2 


(28) 


Note that Mi = (nib + mt)I and M 2 = [**'^** are 
constant, and Mq contains the critical cross-compartment 
interaction, by way of which we can use our tail actuator 
(formally acting on an attitude DOF, (f) 2 ) for energizing 
the shank DOF, 62 . 

Let the forward kinematics of the leg be g : 0 ^ 

The constraint in the stance contact mode is 


h ^ : q 


(01-I-(/)i ,02,a;, 2 )^ 



(31) 


and observe that h ^(q) = (s, a) is reminiscent of SLIP 
(§3) and attitude (§4) coordinates. Define 


7rs:=[/4 0]h \ 7ra:=[o/2]h ^ (32) 


The equations of motion are generated in the new coor¬ 
dinates, 

y = H-^M^(T - N) - H-i(M^C + A^^A)Hy. 

(33) 


In stance, 


ai(q) = ] - R{(j)i)g(9), (29) 

'^Even though the dynamic task here is quite different from free- 
fall, in the language of [13] this is saying that the tail should be light 
but effective. 
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where ^ := 9i — (j )2 (the tail-leg angle), and the right 
summand in (34) is quite clearly the disturbance caused 
due to the added attitude degrees of freedom. 

With the same choice of H, we can similarly recover 
weakly decoupled flight dynamics: 




sin(0i-|-(/)2) 

— COs((/ll-|- 02 ) J ’ 


(36) 

(37) 


5.4 “Physical” Decoupling and Anchoring 

With the highly restrictive assumption 5 (allowing for 
infinite tail inertia), the tail motion is essentially negligi¬ 
ble. Under these conditions, we show the emergence of 
the beginnings of a classical anchoring relation [ 10 ], via 
a natural (weak) decoupling of the 6 DOF dynamics into 
“point-mass” and attitude compartments. A more general 
analysis that is more physically relevant is forthcoming in 
future work. 


Proposition 8 (Flow-invariant submanifold). Under as¬ 
sumption 5, in each hybrid mode, (i) the submanifold 
'Ll = {Tq e TQ ; Tfi = T(j )2 = 0} is invariant un¬ 
der the action of the flow generated by f!["^, and (ii) in 
each hybrid mode, the closed-loop flow restricted to U, 
Tq = fl'^{Tq\ii) is a cross-product of the template vec¬ 
tor fields, 

fr = ffo^sX f^OTTa, (38) 


where tTs cind tTq represent projections to the SLIP and 
attitude components o/q respectively. 


Proof Applying assumption 5.ii to the equations of mo¬ 
tion, the plant dynamics Tt)) are 


I stance 


Tf, 20263 



1 Tt [ sin C/02 

ks(Pl-e2) g2 

L " J 

ptruf, [ -cosC 


distance [ ] ? 

[f ] Iflight = [-^g] + 
fl|flight = [rt ] ) 


sin(0i-|-<?i2) 

-COs(</)l-|-(/)2) 


(39) 


We can check that we have available affordances 
through our two actuators to assign (scaled versions 
of) our template controllers in Table 2, (i) r/i| stance = 
—( 7 ]’(ai,di) to control ai, and r/i|flight = ^ con¬ 

trol X, and (ii) rt|flight = 52 *^( 02 , 02 ) to control 02 , and 
distance = -ptd 2 mb • 5 i(i) to contTol hopping height^"^. 

Under assumptions 3.iv and 5.iv, we show that the high¬ 
lighted terms in (39) vanish inside U: 

'^We observe that by assumption 3.ii, 6^2 pi is roughly constant, 
so the scaling need not be configuration dependent. 


i) M 2 ^ 00, so in the dynamics equations d = 0. Re¬ 
stricted to U, a = 0. This proves part (i) of the claim. 

ii) From d = 0 and (21), r^l stance = A | flight = 0. 

iii) Since (?)2 = 0, ^ ^ 0 (from assumption 3.iv). 

By comparing the thus-restricted plant dynamics (39) 
to (15), (16) and (20), we obtain part (ii) of the result. 

□ 

Additionally, the invariant submanifold in the flow 
leads to an invariant submanifold in the hybrid execution: 

Proposition 9 (Return map-invariant submanifold). The 
set li is invariant under the return map T'''^(Tq|u), and 
restricted to U, o tts x o 

Proof We first define the return map F‘^'^ by instantiat¬ 
ing a “cross-product” hybrid system U"^) as 

(a) 2 ?^'" := T" x (b) := F x P, and (c) as 

defined in Proposition 8 , where Tf := TS^ x for each 
i (ensuring = 0) and r? : T? —> is defined 

^ [i7rmod27r] ' (^0) 

With these modifications, the fa dynamics (22) are ig¬ 
nored, and the clock of the HIR subsystem is being driven 
by the SLIP subsystem*^. This ensures that the conditions 
of Proposition 7 still hold, i.e. vTa o F^"' = F^ o TVa- 
Additionally, the decoupled nature of (Proposi¬ 

tion 8 ) allows us to conclude that tt^ o F^"' = F^ o tTs, so 
that 

= TT, O XTTaO F^"' = T" O 7r« X T" O TTa, 

which concludes the proof. □ 

We leave to future work a proof that U is attract¬ 
ing, which is a requirement for demonstration of anchor¬ 
ing [ 10 ]. 

6 Experimental Results 

In this section we present empirical data obtained from 
the Jerboa (§5.1). In the first three subsections, we present 
data from a few “nodes” of our composition tree (Fig. 1). 
Finally, a crucial examination of our idea of composition 
of templates, when implemented on the Jerboa, is pre¬ 
sented in §6.4. 

**This coupling interaction importantly invalidates the uia- 
dependent bound on k (26). Our solution is to scale the input such 
that k is high enough for the shortest feasible transition time in verti¬ 
cal hopping. 


















Figure 12: Two datasets corresponding to different tail masses: The blue traces use the rrit = 150 g (as in Table 4), but the red 
traces use mt = 100 g. Note that the tail displacement is larger for the lighter tail mass, although vertical behavior is largely 
unaffected. 


6.1 Effect of Varying Tail Mass on Vertical Hop¬ 
ping 

The first empirical result we present corresponds to the 
top left leaf of Fig. 1 —empirical vertical hopping. In or¬ 
der to facilitate the analysis in this paper, in assumption 5 
we stipulated an ideally effective [13] tail, with negligible 
mass and infinite inertia. We connected the robot (Fig. 7) 
to a boom and constrained the body pitch as well as the 
fore-aft DOF. By varying the tail mass (with a fixed tail 
length given in Table 4), we obtained two vertical hop¬ 
ping datasets plotted in Fig. 12. 

We observe the following: 

i) Increasing tail mass results in smaller tail displace¬ 
ments. Taken to the limit, this sheds some light on 
assumption 5 : a large tail mass would indeed render 
the tail motion negligible. 

ii) The hopping height remains relatively unchanged in 
spite of this physical variation. From (34), the force 
acting on the leg-spring depends only on the (feed¬ 
forward) tail torque, Tt (as in Table 2). 

Consequently, we see that the tail mass is a tunable de¬ 
sign parameter that allows us to trade off the conditions 
of assumption 5 (negligible mass versus large inertia— 
both affecting coupling interactions) without affecting the 
vertical behavior. 

6.2 Empirical Validation of Attitude- 
Decoupling Change of Coordinates 

An important foundation of our attitude control strat¬ 
egy is the decoupling of the two attitude DOFs (§4), such 
that ai is controlled in stance, and 02 is controlled in 
flight (20). However, the body pitch and tail angle are 
clearly coupled in flight*^. To resolve this, as shown 

'^Since the tail actuator is attached between the body and the tail, 
tail torques are felt by the body. 


in (31), we use as a decoupling change of coordi¬ 
nates. 

In terms of implementation this strategy requires the 
estimation of a single scalar parameter that defines M 2 up 
to scale (see the text just after (28)). To test our the change 
of coordinates empirically, we suspended the robot about 
the CoM and applied a feedforward sinusoidal Tt signal. 
The resulting traces for the physical attitude coordinates 
are shown in Fig. 13. 

Recall from (37) that in flight, di = 0. In practice, 
we observe from the right of Fig. 13 that there are small 
ai-variations are at a much slower time scale than 02 - 
variations. The reason that di is not zero is that we were 
unable to suspend the robot at precisely the CoM, and so 
gravity exerts a net moment on the body—appearing as 
a slow ai-oscillation. Other than this minor deviation of 
our physical platform from assumption 5, it appears as 
if the attitude-decoupling change of coordinates is indeed 
effective. 

6.3 Trading off Forward Speed and Hopping 
Height for “Leaping” 

The “stepping” fore-aft control using the touchdown 
angle as a control input (10) essentially allows us to trade 
off vertical and fore-aft energy—appearing as a pure ro¬ 
tation in (12). Even though for steady-state behavior we 
choose the touchdown angle to stabilize forward speed, 
it also allows for transient behaviors such as a one-shot 
“leaping” motion (term coined by Raibert [14]). In partic¬ 
ular, choosing a larger (in magnitude from vertical) touch¬ 
down angle than that dictated by (10) results in added ver¬ 
tical height and reduced fore-aft speed. 

The results of an empirical test of this one-shot leap¬ 
ing strategy are showing in Fig. 14: we can indeed get 
a large increase in apex height using this strategy. This 
kind of “asymmetry” [14] or deviation from steady-state 
may have important applications in behaviors that require 
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Figure 13: Testing our decoupling change of coordinates from the physical body pitch, tail angle coordinates (^i, (j> 2 ) to our 
chosen attitude coordinates (oi, 02 ) by suspending the robot about its CoM (see §6.2). 


rapid changes in the body energy, and we plan to explore 
more such behaviors in future work. 

6.4 Empirical Validation of Composition 

By physically constraining some of the DOFs, we test 
our hierarchical composition (Fig. 1 ) at as many “nodes” 
of the composition tree as possible. Note that it is infea¬ 
sible to isolate the fore-aft or the closed-loop pitch cor¬ 
rection templates in a physical setting. The results are 
summarized in Fig. 15. Five strides are averaged within 
each category, and aligned with ground truth knowledge 
of the touchdown event. We observe that 

i) there is a vertical limit cycle that retains its rough 
profile and magnitude through three anchoring bod¬ 
ies, 

ii) the hip angle roughly satisfies = 0 in sfance and 
fhe sfance duration is roughly consfanf (corroborafing 
assumptions 3.ii-iii, and our MBHop model (11), 

hi) fhe shape coordinafe is desfabilized in sfance and sfa- 
bilized in flighf, and fhe pifch-deflecfions are small 
in magnifude over fhe sfride, and in agreemenf wifh 
( 22 ). 

Qualifafively, fhe “failed poinf-mass hopper” configura¬ 
tion affained sfable forward hopping af confrolled speeds 
upwards of 20 sfrides, only limifed by space. The fully 
unlocked sysfem has so far hopped for abouf 10 sfrides af 
mulfiple insfances before failing due fo accumulated error 
causing large deviations from fhe limif cycle. We believe 
fhe prime reason for fhis is fhaf fhe CoM is significanlly 
afl of fhe hip (violating assumption 5.i). We aflempfed fo 
compensafe for fhis effecf wifh a counferbalance visible in 
Fig. 7, buf an unaccepfably large weighl would have been 
required fo complefely correcf fhe problem. 

In fhe video affachmenf, we include clips of fhe robof 
hopping along a boom, wifh varying degrees of physical 
consfrainf corresponding fo fhe “bodies” of Fig. 15 (an- 
nofafed in fhe video). The confroller implemenfed on fhe 


hardware is agnosfic of fhe physical consfrainf, and fakes 
fhe decoupled form of a cross-producf of fhe rows of Ta¬ 
ble 2. 

7 Discussion and Conclusion 

Raiberf’s hopper [14] made significanf empirical ad¬ 
vances in fhe field of robofics, buf fo our knowledge, no 
previous accounf in fhe liferafure has provided any formal 
condifions under which such simple and decoupled con- 
frol sfrafegies will work. In fhis paper, we apply simple 
decoupled confrollers using similar ideas (including fhe 
exacf same fore-afl (10) and pilch (22) confrollers), buf 
wifh a new vertical hopping scheme (§3.1) and a new fail 
appendage fo enable if. Moreover, we consfrucf absfracf 
models (fhaf appear fo, neverlheless, be represenfalive of 
empirical dafa) fhaf enable us fo presenl analyses of slabil- 
ify for each of fhese subsystems, and make sfeps fowards 
a local proof of sfabilify for fhe failed hopper (a subjecf of 
fulure work by fhe aufhors). 

The firsf focus of fulure work is a complefe analysis 
of sfabilify of fail-energized hopping on fhe Jerboa, and 
developmenl of formal fools for design and verification of 
parallel composilion. Second, our analysis in Ibis paper is 
very specifically fargefled fo fhe failed hopper (including 
fhe hand-designed hierarchy in Fig. 1), buf in fulure work 
we plan fo generalize fhese ideas fo olher lasks as well 
as plalforms. As explained in §2, we focus on closed- 
loop lemplales in fhis paper, buf Ihere is an accompanying 
inferesling problem of assignmenl of aclualor affordances 
fo fhe conlrol of specific compartmenls. ar 
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Figure 14: Top: Snapshots of fore-aft hopping behavior in a trial where we test a “leaping” motion [14] —the robot stubs its toe 
at the last touchdown in order to gain a boost in vertical height at the expense of forward speed (see §6.3). The red line shows the 
CoM-trajectory of the robot. Bottom: Corresponding traces showing near-steady-state behavior in the fore-aft compartment (leg 
angle, Qi and vertical height, z are plotted) before the “stubbing” event (red overlay). The leg angle shows the “neutral angle” 
with a thin horizontal line, and in order to leap; note that a much larger (in magnitude from vertical) touchdown angle is chosen 
in order to leap. The leg height (z) plot shows the robot getting around 50% larger apex height in the subsequent flight phase. 
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